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ABSTRACT 

The  regression  model  Y  =  XB  +  e,  with  e**»N(0,  o*I ) ,  has 
been  studied  extensively.   That  is,  the  model  in  which  the 
errors  are  independent  and  identically  distributed  as 
N(  O }  fr*)  has  been  studied  already. 

In  this  thesis  we  study  the  model  in  which  the  sample 
observations  are  correlated  with  a  prescribed  correlation 
structure  and  show  that  many  of  the  results  available  for 
the  independent  case  apply  equally  well  for  the  correlated 
samples . 

We  shall  find  that  some  results  obtained  here  are  not 
just  the  same  as  the  case  where  the  errors  are  independent 
and  identically  distributed  as  N(  O,^1). 


TABLE  OF  CONTENTS 


I.   INTRODUCTION 6 

II.   SUMMARY  OF  KNOWN  RESULTS 9 

A.  SOME  PROPERTIES  OF  £  AND  £•*  UNDER  NORMAL         9 
THEORY 

B.  SOME  THEOREMS  USED  IN  DERIVING  THE  RESULTS 10 

III.   SOME  PROPERTIES  OF  B  AND  Y'CY 

UNDER  PRESENT  MODEL 12 

A.  EXPECTED  VALUE  OF  B 12 

B.  DISTRIBUTION  OF  Y'CY^ 12 

C.  EXPECTED  VALUE  OF  Y'CY 13 

D.  DEPENDENT  RELATIONSHIP  BETWEEN  B  AND  Y'CY 13 

IV.   SOME  PROPERTIES  OF  R.  AND  Y   C  Y   UNDER  THE 

NEW  MODEL 15 

A.  EXPECTED  VALUE  OF  B± 16 

B.  DISTRIBUTION  OF  I  £  1/^ 16 

C.  EXPECTED  VALUE  OF  Y*'c*Y* 16 

D.  INDEPENDENT  RELATIONSHIP  BETWEEN  B,  ATTO 

Y  C  Y 17 

V.   INTERVAL  ESTIMATION  AND  PREDICTION  INTERVAL 18 

A.  COVARIANCE  OF  B- 18 

B.  INTERVAL  ESTIMATION  OF  r'B. 19 

C.  PREDICTION  INTERVAL  FOR  ONE  FUTURE 
OBSERVATION 19 

LIST  OF  REFERENCES 23 

INITIAL  DISTRIBUTION  LIST 2k 


I.   INTRODUCTION 

Regression  analysis  is  designed  for  situations  where  a 
variable  is  thought  to  be  related  to  one  or  more  other 
measurements  made,  usually,  on  the  same  objects.   A  pur- 
pose of  the  analysis  is  to  use  data  (observed  values  of 
the  variables)  to  estimate  the  form  of  this  relationship. 

The  problem  we  consider  is  that  of  estimating  or  pre- 
dicting the  value  of  a  dependent  random  variable  y  on  the 
basis  of  some  known  measurements  of  an  independent  controlled 
variable  x.   Scientists,  economists,  psychologists,  and 
sociologists  have  always  been  concerned  with  the  problems 
of  prediction.   Meterologists  are  constantly  analyzing 
data  in  hopes  of  predicting  or  forecasting  with  a  high 
degree  of  accuracy. 

An  example  would  be  to  use  information  on  weight  and 
height  of  people  to  estimate  the  extent  to  which  a  man's 
weight  is  related  to  his  height.   If  among  people  one  was 
picked  at  random  we  might  expect  his  weight  to  be  a  +  bx 
(where  x  is  height),  namely,  E(y)  =  a  +  bx  (where  y  is 
weight).   In  gathering  data,  the  weight  of  every  man  with 
height  x  will  not  be  exactly  a  +  bx.   Therefore  the  differ- 
ence, y±   -   E(y^)  ,  can  be  written  as 

ei  a  yi  -  E(yi)  =  yi  -  a  -  bx«  , 

hence  yi  =  a  +  bxj_  +  ei  . 


The  man's  age  can  also  be  considered  to  be  a  factor  affecting 
weight.   Then  the  model  is  extended  to  be 
•  E(y)  =  a  +  b<\x<\    +  b2x2  . 

More  generally,  we  can  extend  the  equation  to  be 
E(y)  ss  a  +  b^x-i  +  b2x2  +  •  •  •+  bkxk 
yi   =  a  +  b^x-  +  b2x2  + 

A  frequently  used  notation  is 

b„x.~  + 


+  bkxk  +  e. 


*i  =  ViO  +  b1Xi1  +  - 


+  b,  x.,  +  e.   , 

k  lk    1   ' 


for  i  =  1,2,-"  ,n,  with  x>f=  1   for  all  i. 
Now  define  the  following  matrices  and  vectors: 


X  = 


x10  X11 
X20  X21 


x  n   x  . 
nO  n1 


1k 


'2k 


nk 


1  x 


1  x 


11 


21 


1  x 


n1 


L1k 


l2k 


nk 


B  = 


0 


'kj 


Y  = 


n 


e  = 


n  ) 


Then  we  write 

Y  =  XF  +  e  ,  with  E(Y)  =  XB    . 
The  model  Y=XB  +  e   ,  with  e~N(0,crl),  has  been  studied 
extensively.   That  is,  the  model  in  which  the  errors  are 
independent  and  identically  distributed  as  N(  0  7  (T  )  has 
been  studied  already. 


In  this  thesis  we  study  the  model  in  which  the  sample 
observations  are  correlated  with  a  prescribed  correlation 
structure  and  show  that  many  of  the  results  available  for 
the  independent  case  apply  equally  well  for  the  correlated 
samples. 

In  this  thesis,  we  shall  examine  the  model 


Y  = 

V  = 

A  = 


XB  +  e   ,  with  e~N(0,V)  ,  where 
*g(A  +  A»)+°t(I-E), 


a1  ai 
a2a2 


a   a 

„  n  n 


a. 


n  ) 


A'  is  the  transpose  of  _A,  a.  (  i  =  1,2,-*-,n)  and  o(  are  posi- 
tive constants,  _I_  is  an  n  x  n  identity  matrix,  and  E_  is  an 
n  x  n  matrix  all  of  whose  elements  are  unity. 

We  shall  find  that  some  results  obtained  here  are  not 
just  the  same  as  the  case  where  the  errors  are  independent 
and  identically  distributed  as  N(C;o~  )8 
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II.   SUMMARY  OF  KNOWN  RESULTS 


A.   SOME  PROPERTIES  OF  B  AND  O"2  UNDER  NORMAL  THEORY 

The  model  Y  =  XB  +  e,  e~N(0,  a*2,  I ) ,  where  the  components 
of  B  and  <j-2  are  unknown,  and  X  is  a  known  matrix  of  rank 
k  +  1.   Since  the  vector  of  errors  e  is  normally  distributed, 
the  maximum-likelihood  method  will  be  used  to  estimate  B 
and  <r  .   Solving  the  likelihood  equation  we  get 


B 


'0 


A 

b 


k> 


=  S'1X'Y  ,   where  S"1=  (  X'X  )"1  , 


and 


A 


-1 


(  Y  -  XB  )'  (  Y  -  XB  )     Y»  (  I  -  Xj3  'X«  )  Y 


n  -  (  k  +  1  ) 


n  -  (  k  +  1  ) 


a      a  2 

Then,  B  and  (3s*  have .  the  following  properties: 

1.  consistent 

2.  efficient 

3.  unbiased 

4.  sufficient 

5-  B  -  N(  B  ,  Ch'l"1  ) 

6.  complete 

7.  minimum  variance  unbiased 

8.  (  n  -  k  -  1  )  <^2 


O- 


% :  (  n  -  k  -  1  ) 


9o   B  and  <?-  are  independent. 


B.   SOME  THEOREMS  USED  IN  DERIVING  THE  RESULTS 
Theorem  1,   If  Y  is  distributed  as  NCM^i),  then  JL'A 


O^ 


2.  f 

is  distributed  as  X(k,  A)  ,  where  ~\-  ^      ,   and  k   rank 
of  A,  if  and  only  if  A  is  idempotent 


">  0^ 


Theorem  2.   If  Y  is  distributed  as  N(X^V),  then  Y'BY  is 
distributed  as  9C  (K,X)  »  where  A  =  —  Al'B^i  and  k  =  rank  of 
B,  if  and  only  if  BV  is  idempotent. 

Theorem  3.   If  Y~N(.U,y),  then  Y'AY  and  Y'BY  are  independent 
if  and  only  if  AVB  =  0. 

Theorem  j£.   If  Y~N(U3\/),  then  C'Y  and  Y'AY  are  independent 
if  and  only  if  C'VA  =  0. 
Theorem  5-   (Hogg  and  Craig  theorem) 

Let  Q  =  Q1+Q2+---  +Qk,  where  Q,QV  •••  ,Qk-1and  Qk  are  k  +  1 
random  variables  that  are  quadratic  forms  in  the  observa- 
tions of  a  random  sample  of  size  n  from  NC/U,^-). 
LetiL  be  %Z(r)  »  let -4  be  J^ir.)  »  i=  1>2'  '  "  »k~1  »  and 
let  Q,  be  non-negative.   Then  the  random  variables  CL  ,  Q?, 
**',  Q,  are  mutually  stochastically  independent  and,  hence, 


<*i 


k-1 


i-is  X^rk=r-?="  ri} 


Theorem  6.   (Baldessari  theorem) 

Let  Y  be  N(i},  V  )  and  B  ,  B  ,  *  **,  B,  be  (  n  x  n)  idempotent 

matrices  satisfying 

2_  B  .  =  I  -  —  E     ,  where  X  an<3  E  are  matrices  as 

defined  above. 

Let  >X  be  a  positive  constant. 
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Then  a  necessary  and  sufficient  condition  for  X'B.Y,   , 
j  =  1,2,'  •  -  ,k  ,  to  be  mutually  independent  and  have 

non-central  Chi-square  distributions  with  r .  (r  .=rank  of  B. 

J    J  J  * 

j  =  1,2,-  •  ■  ,k)     degrees  of  freedom  is  that  the  covariance 
matrix  V  has  the  following  structure 

V  =  5g(A  +  A'  )  +  d  (  I  -  E  )  . 
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III.   SOME  PROPERTIES  OF  B  AND  Y'd-XS  1X')Y 
UNDER  THE  PRESENT  MODEL 


We  consider  the  model 

Y  =  XB  +  e  ,   e~N(0,V),  where  V=^(A+A')+*  (I-E)  . 
Let   B=S"1X»Y   and   X'CY    (Y-XB) » (Y-XB)     Y" (I-XS*1 X« )Y 


Now  we  will  show  that 

A.  E(B)  =  B 

B.  Y'CY     _2 

5Tgasa  —  X  (  n  "  k  ~  1  ) 

C.  E(Y»CY)  =  c<  (h-k-1  )   and 

D.  B  and  Y'CY  are  not  independent. 

A.   EXPECTED  VALUE  OF  B 

E(B)  =  E(S"1X'Y)  =  S"1X'E(Y)  =  S~1X'E(XB  +  e  ) 

=  B 
Hence  B  =  S  X'Y  is  an  unbiased  estimate  of  B 


B.   DISTRIBUTION  OF 


(  Y  -  XB  ) '  (  Y  -  XB  ) 


We  have  seen  that  (n-k-1  )  >"  /yz(n.k.1  )   f  where  fc\  _iLSJ. 


gs.-  n-K- 1 

-1 


Therefore  we  can  conjecture  that  —  v—  —  —  -   is  distrib- 

uted  as  )(_"(n-k-1  )  .   The  necessary  and  sufficient  condition 
for  X'CY   to  be  Xz(n-k-1  )   is  that  2  CV  is  idempotent. 


*  * 


Now 


we  want  to  show  that—  CV  =  —  (I.-XS  2£')X.is  idempotent, 


that  is,  to  show  that  ±(I-XS~1  X'  )V(I-XS"1  X'  )  =  I-XS~1X»- 


-^ 
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Note    -i-(I-XS-1X')V(I-XS-1X') 

=i  (i-2^"1x,)(1^(a+a')+  (I-E))  (I-XS~1XM 
=  (i-xs"1x»)  , 

since  X_S"1X*E=E  ,  XS~1X»A'=A»      and  AXS~1X'=A  . 
Hence,  by  theorem  2,   X'2X  ^.^O  q  ,  A  ) 

where  q  =  rank  (!-%$■     £f  ) 

=  rank  (  I  -  XS  X'  ) 

=  tr  (  I  -  X_S"1X'  ) 

=  tr  (  I  )  -  tr  (  X'Xjf 1  ) 

=  n  -  k  -  1 

since  I-XS"  X*  is  itself  symmetric  and  idempotent, 

_-] 

and   X  =  kptyt(  -  Z   —     -'  )  XJ3  =  0  . 


Therefore    Yf CY  ^  yZ   j  n_k_i  ) 


C.  EXPECTED  VALUE  OF  Y'CY 

From  B  it  is  obtained  directly  that   E(Y» (I-Xjf  X1 )Y)= 

c*(n-k-1  )  also  then  -^-^     is  an  unbiased  estimate  of  oC  . 

n-k-1 

D.  B   =   S_1X'Y   AND   Y'CY   ARE   NOT   INDEPENDENT 

To   show  that   they  are   not   independent,    let's   show 
S"1X'VC   =   S"1X/VjI^XS~1X* )    ^   0    (    Th.    i|    )    • 
S"1X'V(I-XS"1X')-S"1X'  (lg(A+A«)+  --UI-E))   (I'M"  V) 

=^(S"1X/A/-S~1X/A/XS"1X»)  *  0  • 

A 

We  recognize  that  B  and  Y'CY  are  not  independent,  which 
is  not  the  same  results  as  obtained  in  the  i.i.d.  case. 
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/N 


For  reference,  the  variance  and  covariance  of  B  is 
provided  as  follows: 


V(B)=V 


/  — 

a 


-1- 


-*(n-1 )/n+(nax,-a'X1+ix»)S   x 
Sj1  (^'a-^nax-ix) 


T\ 


where  a  =  (1/n)E  a.  ,   a'  =  (  a1  ap 


X!  =  (  X..   X. 


(^a»X1-^nax,-«xx»)S1 


-  a   )  , 
n   ' 


-1 


1  x,2  '"*  A#k 


),  x.  .=  (1/n)5ixi  .  for  all  j  , 


i  =  i 


Si- 


and 


'11 


n1 


"1k 


nk 


S-  =  XJ^  -  nxx' 


Note  that  c<  factors  from  the  lower  right-hand  corner, 
"but  not  from  the  rest  of  the  matrix.   This  causes  trouble 
"because  the  unknown  parameter  cA  does  not  cancel  out  in 
forming  the  usual  Student  t  statistic. 

Therefore  we  develop  the  new  modified  model  without 
Bq  in  the  next  chapter. 


U* 


IV.   SOME  PROPERTIES  OF  j^  AND  Y*  (  I  -  X*S*"1X*  )Y* 

UNDER  THE  NEW  MODEL 


We  consider  the  model 

Y  =  XV"^1  +  e   , 
where   Y*  =  Y  -  y1_  =  (  I  -  (1/n)E  )Y 

x*  =  X-i-Cx.-ji:  '  '  •  :*'kl)  =  (  I  -  (Vn)E  )XJ] 

B^    s  (  b1  b2  •  .  .  bk  )• 

e*  =  e  -  51  =  (  I  -  (1/n)E  )e   . 


Note  then  that  y.  =  y.-y  ,  x,  .  =x,  .  -x.  .  and  thus 

£  y*=  iz  xtV0'  that  is>  1  ,Y*=0'  1'^=  °  • 

Then 

V(e*)  =  oL  (  I  -  (1/n)E  )  , 

i.e.,   e*~  N(  0  ,  oi(I-(1/n)E  )  ) 

and      Y*  ~  N(  X'%  ,   -A  (  I  -  (1/n)E  )  ) 

under  our  original  assumptions. 

Let  8.=  S   X*  X*"   '  where  S"  =  X"  X"  ,  and 

Y*V*Y*  _   Y"^I-X'V:'"1X""*')Y*""  =   (Y'""-X'""B1  )  '  (Y^-X*^  ) 


We  are  going  to  shov/ 

A.  The  expected  value  of  E^ 

B.  The  distribution  of  I  £  I 

•w-  '  -5C-  -£ 

C.  The  expected  value  of  X  C  Y 

D.  That  B1  and  Y  C    ?      are  independent. 
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A.   EXPECTED  VALUE  OF  B^ 

E^)  =  E(S"  nX  Y*) 

-"--1  -"-»    -* 
=  S"  V  E(Y  ) 


=  S*"1X*  E(X*Bt+e*)  =  a, 


A.      -"--1   #'   "& 

Hence  B.=S   X  Y     is  an  unbiased  estimate  of  B.  it  is 
easily  shown  that  B.  is  the  same  as  B.  in   £  ~  >su 

B.   DISTRIBUTION  OF  (VoO(Y  C  Y  ) 

Since  £*"\^  (I-is)=I-X'""s'""1  X*"'-^  E    is  idempotent, 
(1/0c)Y'"',C*V"^(q,X)  by  Th.  2,  where 
q  =  rank  (lA)(I-fs*"1X*') 
=  rank  (I-X"S*"'"1  X*'  ) 
=  tr  (I-X"S""nr  ) 
=  tr  (I)  -  tr  (X*fX*S*"1) 
=  n  -  k 
and   X  =  ^B1'r^1/^)(I-xV'"1X"^X%  =  0 


Therefore 


(iAt)(Y*fcV")  ^%  ( ^  - k ) 


t 


C.   EXPECTED  VALUE  OF   Y"  C"Y" 

It  is  obtained  directly  from  B  that  E(Y  C"y" )=  oC  (n-k) . 

-"- '  -'-  -::- 

•         •              v   c  Y 
Also  it  is  easily  seen  that  = ~  —        is  an  unbiased  estimate 

of  "C . 


16 


D«   Sis:5*""12*,I*    AND  I^'c^Y*   ARE  INDEPENDENT 
S^^X^'ocCl-d/njEjCl-X^S*"1^1) 
=  o((S   X  -(1/n)S""nX"  E)(I-X  S   X   ) 
=  0        (  note  X*'e  =  0  )   . 

Therefore,  by  Th.  k,    they  are  independent. 


1? 


V.   INTERVAL  ESTIMATION  AND  PREDICTION  INTERVAL 

The  construction  of  interval  estimators  and  a  prediction 
interval  to  contain  a  single  future  observation  in  the  i.i.d. 
case  with  a  regression  model  has  been  discussed  in  many 
standard  texts.   In  this  chapter  we  extend  these  results 
to  the  case  where  the  samples  are  correlated  with  a  special 
type  of  covariance  structure.   We  shall  discuss 

A.   Covariance  of  B,  ,   V(B,.  ) 


B.  Interval  estimation  of  r'B, 

C.  Prediction  interval  for  one  future  observation. 


S\  A 


A.   COVARIANCE  OF  B      ,   V(B,  ) 

Since  B.  is  equal  to  the  product  of  a  constant  matrix 

#--1  *•»  .  * 

S   X    and  a  normally  distributed  vector _T* ,    B.  has  the 

k-variate  normal  distribution. 

We  have  already  shown  that  the  mean  of  B.  is  B« . 

a 

The  covariance  matrix  of  B1  is 

V(|1)=E((B1-B1)(|1-B1)') 

=E(  (S*"1X*' £*..&,  )  (S*~1X* 'y*-^  )  '  ) 

=e(  (s*"1 x*'  (x%  +**)-£,  )  (s*-1  x*f  (X%  +£*)-&,  )  '  ) 
=S   'X  V(e  )X  S 

=  *  S'""1x'"'f  (I-(1/n)E)X*S*~1 

-"--1 

So  B«  is  distributed  N(  B.  ,  k   S*   ) 
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B.   INTERVAL  ESTIMATION  OF  r'B^ 

Frequently,  an  experimenter  is  interested  in  setting 
confidence  limits  on  some  function  of  b,  .   A  method  is 
available  for  one  of  the  most  frequently  occurring  cases, 
i.e.,  for  the  case  of  a  linear  function  of  the  b.  . 

We,  here,  discuss  the  case  of  a  linear  function  of  the  b. 

(except  b0  ).   Let  r  be  a  known  k  x  1  vector  of  constants, 

then,  to  set  confidence  limits  of  size  l-d\,    on  rB.  ,we 

proceed  as  follows. 

Since  r'B.   is  distributed  N(r 'B, ,xr ' S"   r   ) ,  it  is 
/«. 
clear  that   g'J^-r'B^      is  distributed  N(0,1).   We  also 

know  Y  C"Y"  ^  'V~(n-k)    and  independent  of  B.,  . 


*. 


h 


Therefore,   I'^-r'B    /  (n-k) 

•7- — Lsn — —  k>  «  *l  ~  t  (  n  -  k  )  • 

Ur'S;   r)^  \y"    C**Y  / 

We  arrive  at  the  probability  equation 

PCr'B.-ctct,  <   r'B.  <■      r'^+cU  )  =  1-  <*j  , 

where  c  =  f    (r  '  S*~1  r  )Y*'  C"'V: 


\     n-k 
C.   PREDICTION  INTERVAL  FOR  ONE  FUTURE  OBSERVATION 

Hahn  /~^J7   has  derived  prediction  intervals  to  contain 
some  future  observations  of  ia.d.  random  variables  from  a 
normal  distribution  for  a  regression  model.   In  this  section 
we  extend  this  result  to  the  case  where  the  samples  are 
correlated  and  have  a  special  type  of  covariance  structure 
to  construct  prediction  interval  to  contain  a  single  future 
observation. 
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Let  y  +1  "be  the  one  future  sample. 

Define  Y  to  be  the  (n+1)  component  vector  including  the 
original,  sample  of  size  n  and  y  +«  •   We  assume  Y~ N(XB,  V) , 
where  X  is  (n+l)x(k+l),  of  rank  k+1  and  V  is  as  defined  in 
III,  and  is  (n+l)x(n+l). 

Define  y*+:.  as  follows: 

yn+1=yn+1^n,1   '  Where   yn+1=  (V(n+1))  gy.   . 
Then,  y*+1=yn+1-(1/(n+D)(yn+1+nyn)  =  (n/(n+1  ) )  (yn+1  -yj  , 

•  where  y  =  (1/n) >_  y. 

n      i=1  x  ' 

and  y*,-,  can  "be  expressed 

y*+1=(n/(n+1))(blX*+1}1+b2x*+1j2+  ...    +bkx^+1)  k+e^+1) 
=(n/(n+1))(   B{X*+1+  e*+1    )      , 

where     x*+1  =    (   x*+1    1    ,    x*+1    2   ,  •-.  ,    x*+1    k   )■    . 

We   also,    therefore,    know  E(yn+1  )=nAn+M  B^X^+l         and 

V(y*+1)=  (n/(n+1))<      since 
V(y*+1)=(n/(n+1))2(V(yn+1)+V(yn)-2C0V(yn+1,yn)) 

=  (n/(n+1))2(an+l+a-((n-1)/n)x-2(^a+^an+1-A)) 

=  (n/(n+1))*     . 


Next,    let      yn+is<n/(n+1>^iS+i 

then    E(y*+1)  =  (n/(n+1))B1'X^+1  and 

VCy^^-otdi/Cn+l))2^;^*-1^ 
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Now  let  z  =  y*+1-  y*+1   , 

then  E(z)  =  0 

V(z)  =  <>c(n/(n4-1))(1+(n/(n+1))(X^1s"-1X^+1))   . 

Therefore,  z-N(  0  ,  *  (n/n+D)  (1  +  (n/n+1))X*^  S*~1  X*+1  )  )  , 

The  standardized  variate  z'  =  z/(V(z))^    is  normally 

distributed  with  mean  0  and  variance  1.   Furthermore, 
since  (1/a)Y  £"Y"~X  (n~k)   and  i-s  independent  of  z  ' , 

the  ratio 

zl        4   a*   /  n  -  k  . 

T=/?Wy"  =(yn+1"y;+1  \((n/n+1  )+<n/n+1  J^S^Cl  ^Vj 
\  *(n-k)/ 

follows  a  Student  t  distribution  with  n-k  degree  of  freedom. 

■*  •  *  ■* 

Proof  of  independent  relationship  between  Y  C  Y  and  z. 

*•    *   *  *  A* 

Y     C   Y     and   z   =  v.-    -   y     -    can  be   expressed   in  terms   of 
—     —  —  •'n+l        Jrnrl 

Y  =    (ylf     "*,    yn+1)'    as   follows: 

Y*'c*Y*=  Y'[l"(Q{n)-)c"(I-(1/n)E;0)Y  =  Y'D'C*DY 

=  (n/n+1)((-l,---i,1)-r)1S""1rT(I-(1/n)E;0))Y 
=  (n/fn+1))BY 

Then  it  can  be  seen  that  BVD'C  D  =  0,  where  V=~>(A+Af  ) 
+  d(I-E)   is  an  (n+l)x(n+l)  matrix. 

Therefore,  by  Th.  4,  they  are  independent. 
Consequently,  a  two  sided  lOOrfo   prediction  interval  to 
contain  the  future  sample  y  +1  using  the  regression  obtained 
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from  the  n  given  observations  is 

,  ((n/n+1))  +  (n/(n+1))2X*^S*-1X*+1  )(f  cV) 


n  -  k 


Finally,  the  lower  limit,  L,  and  the  upper  limit,  U,  for 
y  +1  are  expressed  as  follows: 

L  =  ((n+1)/n)(y*+1-  t r  u  )  +  yn 
U  =  ((n+1)/n)(y*+1+  t £  u  )  +  yR 

/l(n/fn+1))  +  (n/vn+1))X*'  S*"1  X*   )  (Y* '  C  V")N 

where  u  =  I UZ1 EZJ 

n  -  k 
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